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SUMMARY 

Given a "new" algebraic equation which can be regarded as a variation 
of an "original" algebraic equation whose well-separated roots are known 
initially, the roots of the "new" equation are here derived explicitly in terms 
of Newton's well-known formula. The adaptation of this formula to cases 
where the roots of the "original" equation are not well separated is also 
considered, and numerical examples are included. The procedure can be 
adjusted satisfactorily whatever the arrangement of the roots of the "original" 
equation may be, and is suitable for use in a digital'computer programme. A 
significant saving of machine time relative to an iterative process dealing 
with the "new" equation in isolation is to be expected. 



1. INTRODUCTION 

In a recent investigation undertaken by the 
author concerning the propagation of electromagnetic 
waves through the ionosphere, it was necessary to 
solve a number of quartic equations having complex 
coefficients which varied slowly (and normally 
continuously) with height; horizontally the iono- 
spheric conditions were assumed to be constant. 

A computer programme (in Elliott autocode) 
was devised for solving such a quartic ab initio, 
and for deducing the relevant eigenvectors associa- 
ted with ionospheric propagation at any height, but 
there is a possibility of saving considerable 
machine time if advantage is taken of the fact that 
the roots of the quartic vary continuously with 
height. 

In straightforward cases where the roots of the 
quartic at the initial height Zn are well separated, 
the roots at a neighbouring height Zq - h can be 
deduced explicitly but approximately by an adapta- 
tion of Newton's method of root approximation. 
This adaptation is explained in Section 2, and its 
limitations are discussed. Alternative procedures 
required when the roots are not well separated are 
considered in Section 3. Extension to equations of 
higher degree is briefly considered in Section 4, and 
Section 5 gives examples. There is thus no essen- 
tial originality in this report; its objective is the 
purely practical one of making other people's work 
easier. 



2. THE CASE OF WELL SEPARATED ROOTS 
Suppose that the original equation is: 



F{x) - x"^ + a„x° + a^x'^ + a^x + ao 

^ (x - \) (x - A2) (x - Xg) (x - \0 ^ 

(1) 

and that the roots \, X?. ^ and X^ are well sepa- 
rated, so that none of the differences X,- - \.j need 
be regarded as small. Now suppose that the coef- 
ficients ai are changed to a; + Sai where the 
changes Sa; can be regarded as small. Then the 
resulting changes SXj in the roots are given (to the 
first order) by Newton's formula: 



S\,- -= 



FKX,) 



^(X,-) 
F\ki) 



where 



iagX^ + ha^x^ + ba^x + Sao 



(2) 



(3) 



and FKx) means the derivative of F(x) with respect 
to X. Equation (2) shows that it is the size of 
+^(X,) rather than that of Sa, that matters. Now it 
can be shown that: 



FWi) - Viiki - kj) 



(4) 



where the symbol U indicates the product of all 
factors (X, - Xy) for which / A i. Equation (4) 
indicates that F^X,) will be small if any of the 
other roots kj is close to Xj, and explains the well- 
known fact that (2) is unsatisfactory in such cases. 
But if the numerator ^(X,) is also sufficiently 
small, (2) can still be used. Furthermore, 4^{ki) 
and the original roots Xj are known ab initio, so 
that when the original roots are well separated. 



FK^i) is not small and S\i is of the same order of 
magnitude as <^(^i). If two roots of (1) are close 
together, so that (ki • Xfe) is small but none of the 

other differences {ki . kj) are, then (p(ki) must be 
of order {ki • k^y to permit the use of (2) as it 
stands. If three roots are close together, so that 
{ki • kk) and {ki - ki) are both small quantities of 
the same order of magnitude, say M, then F\ki) is 
of order M^from (4) and therefore 4^{ki) must be of 
order M° to permit the use of (2). 

When Equation (2) can be used, {ki + Ski) is a 
closer approximation to a root of: 



Fix) + 4^{x) - 



(5) 



than ki was, and the procedure discussed above can 
be repeated so that the roots of (5) are obtained to 
any required degree of accuracy. For suppose we 
are seeking a better approximation to the root near 
\ + SAi and that; 



^{x) - 



4Kx) 



{x-k,){x-K)(x-K) 

[so thati//i(\) - •i>{\)/F\K)]. 
Then (5) can be written: 

(X - \) + s^iCx) = 



(6) 



(7) 



and Equation (2) (for i = 1) is equivalent to replac- 
ing \ii(x) in (7) by ^^^{ki). If <;6(x) is of the first 
order of small quantities, then so is 0i(x) in the 
case of well-separated roots. We can therefore 
apply Taylor's theorem and write: 



^i(x) W 0i(\) + (x - \) ^\{\) 



(8) 



Substituting from (8) into (7), we thus obtain: 

X - \ - S\ - X - \ + V>i(\) - -{x - K)^\{K) (9) 

If in (9) X is the root of (5) nearest \, (x - \) is of 
the same order of magnitude as 8\ = -^/'i(\) and 
therefore (x - \ - S\) is of order (S\) ^. Hence in 
the case of well-separated roots, the correction 
required to the approximation given by (2) is an 
order of magnitude smaller than the correction given 
by (2) itself to the initial approximation \. (The 
argument is equally applicable to any other root ki, 
but the verbal expression of it is more complicated.) 

In cases mentioned above where the roots are 
not well separated, but (2) can nevertheless still be 
used, the corresponding result can be shown to be 
that the correction given by (2) is of the same order 
of magnitude as the separations of the roots which 
are close together, while the correction subse- 
quently required is of the order of the squares and 
products of the separations of the roots. 



3. ALTERNATIVE PROCEDURE WHEN 
ROO TS ARE NOT WELL SEPARATED 



THE 



The case in which the foregoing argument is 
unsatisfactory is that in which the roots X.( are not 
all well separated, and where <p{ki), although small 
in the absolute sense, is not sufficiently small 
relative to FKki) to make S\,- small in Equation (2). 
Consider first the simplest case of this kind, when 
Xi and A2 are nearly equal but are well separated 
from Xg and k^. In such a case divide through (5) by 
(x - Ag) (x - \^ and let: 



^Pl2(x) - 



4>{x) 



(X - Ag) (X - \;) 



Then instead of (7) we have: 

(x- \)(x- As) +<^i2(x) - 



(10) 



(11) 



which can be rearranged in the form: 

{x-'/2(Ai + As)}^- { :2— 3)'-^,2(x) (12) 



.[^).,.. 



Now in the case of well separated roots con- 
sidered above, the first term on the right hand side 
of (12) is large compared to the second, whereas 
when A, »j Aj the two terms on the right hand side 
of (12) may be comparable or the term ^i^ix) may 
prevail. In any case we can replace ^'i^ix) by 
i/'i2{''i(\ + Ag)} and then obtain the second approxi- 
mations to the roots of (5) corresponding to A^ and 
Aj as: 

Ml, Mj ^ 
./,(A, . A,)- ± j /^y - ^.i^V'/^ (13) 

If the term i/'i2'f'/2('Ni + Aj)} prevails, the change in 
\ and A2 is of order ['^i2{'^/i(K + k^}'\ '^^ (instead 
of having the same order as 4Aki) as was found in 
the well separated roots case). 

A possible way of getting a better approxi- 
mation to iWi is to replace it by: 



Ail - '/2(Ai + A2) + 



and similarly to replace M^ by: 




A2(Mi)r^ (14) 



'/2 



/^5 = 'MK + K) 




^2)} 



(15) 



and this process has the same sort of convergence 
that is associated with repeated application of 
Newton's formula to the well separated roots case 
as explained in connection with Equation (9). 



Alternatively, having obtained second approxi- 
mations corresponding to (2) for the roots which are 
isolated and to (13) in the case of pairs of nearly 
equal roots, we can regard the equation (5) for which 
a solution is actually required as a variant of the 
equation whose roots are the best approximations 
hitherto obtained, instead of as a variant of the 
original equation (1). If this is done the procedure 
for well separated roots can normally be used there- 
after. If m roots are clustered (m> 2) it may occa- 
sionally be necessary to use the clustered root 
procedure for n roots {n< m) subsequently, as in 
Example 3 below. 

If the quartic equation (1) has two pairs of 
roots such that the members of each pair are close 
together, and each pair is well separated from the 
other pair, we can proceed as in the case of a 
single close pair. The procedure is affected by the 
separation or closeness of roots near the root being 
sought, but not by the separation or closeness of 
roots remote from the root being sought. If the 
quartic (1) has three roots close together and well 
separated from the remaining root iK), (5) should 
be rewritten in the form: 



{x - u)'' - (-1) 



^ Z, say 



(x -K)(x- K) (x - K) 

Now let: 



X - V 



Ag ■= u + r]s 



(17) 



X - u ' y 

A2 - W + ')72 

where 

u - V3(Xi + A2 + Ag) so that t]i + 7^5, + t]^ = 
Then 
(x - \) (x - A2) (^ - \3) 

:= y3 + (7]2T7g + 773771 + T]i772)y - 77^772773 (18) 

If we substitute back into (16) and take all terms 
except y^or (x - u)3over to the right hand side, (16) 
becomes: 



(x - uY 



'^I'^s'^s - ('Hs'^s + ■^s'f]i + '^I'^s) (x - u) - X(x) 

(19) 
ViViVa - X(u) when X "^ u 



(19) thus gives us second approximations to the 
three roots, and the best procedure thereafter is to 
regard the given equation as a variant of the equa- 
tion whose well-separated roots are these second 
approximations,^ instead as a variant of the original 
equation. Clusters of m roots (m> 3) can be treated 
similarly. 



4. EXTENSION TO EQUATIONS OF HIGHER 
DEGREE 

If the original equation corresponding to (1) has 
a higher degree n, but has well separated roots, 
Equation (2) is still valid, (except in so far as <^(x) 
now has n terms instead of four), and so is equation 
(4) except that there are now (n - 1) factors instead 
of 3. 

If the original equation has a cluster of k roots 
\, Xj ... kjc which are close together and have 
average value u and actual values (u + 77^), (u + 77,) 
... (u + 77/c) where ''71 + 772 + ... + 77^ = 0, then the 
second approximation to the roots of this cluster is 
given by: 



k + 1 



"niViVa—Vk 



4iu) 



{u- kk + 1) iu- \k + 2) 



= - X(x), say (16) 



{u - \n) 

(20) 
so that: 
x-u- |z|l/*eKargz ^ 2r-n)/k^^ ^ 0,1 ... fe - 1) 

When second approximations A^, /*4 ^'n to all the 

roots have been found in this way, the given equa- 
tion should be rewritten in the form: 



ix-K)(x- K) ... (x-XA) + n(x) - 



(21) 



and thereafter it can usually be treated as an equa- 
tion having well separated roots. An exception 
occurs in Example 3 below because we have at one 
stage a complex approximation to a root which is in 
fact real. 

In the situation discussed above, we have had 
the advantage that an obvious good initial approxi- 
mation to the roots of the equation to be solved is 
available. When an equation has to be solved ab 
initio, however, the first task is to find such an 
initial approximation. Various well known methods 
are available for this but when a whole system of 
equations with varying coefficients has to be solved, 
one member of the system can usually be found for 
which solution is specially easy. 



5. NUMERICAL EXAMPLES 

In the illustrative numerical examples which 
follow, both the "original" equation and the "new" 
equation with slightly different coefficients whose 



roots are being sought are chosen to have exact 
roots which are stated initially, so that the manner 
in which the approximations derived approach the 
time values can be clearly seen. All numerical 
examples are quartics. 



Example 1 (well separated roots) 
Starting with the quartic equation: 
ix - 0-5) ix - 1) (x - 2) (x - 4) 

-^ x'^- 7-5x3+ 17-5x2- 15x + 4 = (22) 

find the roots of: 

x"- 7-7x + 18-45x2- 16.263X + 4-6494 -- (23) 
of which the exact values are 0-6, 0-9, 2-1 and 4-1. 

In the notation of Equations (1) and (3): 
flg - -7-5, as - 17-5, % - -15, flo ^ 4 
Sog '- -0-2, Sag ^ 0-95, So^ -- -1-263, Soq ^ 0-6494 
\ '- 0-5, Ac = 1, \b - 2, V - 4 (24) 





Equation (2) gives: 




s\ 


0-2304 

- + . 

0-5 X 1-5 X 3-5 


0-2304 
= ^ 0-0878 

2-625 


SX5 


0-1364 


'- -0-0909 


0-5x1x3 


SAg 


0-3234 

- + 

1-5x1x2 


- 0-1078 


sv 


2-0026 

- + : 

"l-S X 3 X 7 


- 0-0954 



(25) 



To make further progress, we regard (23) as a 
viiri ;int ot: 

(x - 0-5878) (x - 0-9091) (x - 2-1078) (x - 4-0954) 

^ X*- 7-7001x^ + 18-452223X'- 16-236464x 

+ 4-612825 :-- (26) 

instead of as a variant of (22) so that the revised 
values of Sa,q etc. become; 

Sag '- + 0-0001 oal = - 0-002223 

Sai ^^ - 0-026536 huo = 0*036575 (27) 



Hence: 



0-020229 



0-020229 



SX^- 



s>i- 



0-3213 X 1-52 X 3-5076 1-713028 
0-01181 

0-010689 0-010689 



0-3213 X 1-1987 X 3-1863 1-227179 



0-00871 



0-028297 



0-028297 



SAi= + 



1-52 X 1-1987 X 1-9876 3-621455 
0-00781 

0-102509 0-102509 



3-5076 X 3-1863 x 1-9876 22-213944 



= + 0-00461 



(28) 



Hence we have all the correct roots with errors only 
in the fourth decimal place after two applications of 
the Newton process. In this illustrative example 
only real roots were involved at any stage. In 
principle no difference occurs when there are com- 
plex roots. For a desk-machine calculation, how- 
ever, the work is more tedious, but this is immaterial 
if a digital computer is used; programming is reason- 
ably straigthforward. 

Example 2 (One pair of roots nearly equal; complex 
roots involved) 

Starting with the quartic equation: 

(x-0-9) (x- 1-1) (x- 2) (x-4) 

= X*- 8x3 + 20-99x2- 21-94X + 7-92 ^ (29) 

find the roots of: 

x'*- 8-1x3 + 22-1025x2- 24-411x + 9-5513 ^ (30) 

of which the exact values are 1*05 + 0-2j, 2-2 and 
3-8. 

liquation (12) becomes; 

0*1728 

(x- 1)2 =0-01 --0*0476 (31) 

3 

so the second approximation for the roots near 1 is 
1 ± 0*218 174j, while the second approximations to 
the other two roots are 2-1714 and 3-8250. Hence- 
forward we therefore regard (30) as a variant of: 

(x2- 2x + 1*0476) (x- 2-1714) (x - 3*8250) 

- X*- 7*9964x3 + 21*3460x2- 22*893029x 
+ 8*700947 = (32) 



and the changes in coefficients are thus reduced to: 

has = -0-1036; Soj ^ 0-7565; 

Sai'- -1-517971; Sa^ = 0-850353; (33) 

SA^ then is found to be 0*0281, which gives the next 
approximation to Xg as 2-1995, which is very close 
to the correct value 2-2. Likewise S\i is found to 
be -0*0237, so that the second approximation to V 
is 3*8013, which again is very close to the correct 
value 3*8. 

To obtain SA^, we can write (30) in the form: 
X- 1 -0-218174J 

-0*1036x3 + 0*7565x^- l*517971x + 0*850353 
(x- 1 + 0-218174J) (x- 2*1714) (x- 3*8250) 

(34) 

On the right hand side of (34) divide the numerator 
by x2- 2x + 1*0476, and replace x by 1 + 0*218174] 
in the remainder and also in the factor x - 1 + 
0*218174] in the denominator. In the remaining 
(real) denominator factors we can either replace x 
by 1 + 0*218174j, or preferably multiply these 
factors out, divide by x^ - 2x + 1*0476 and then 
replace x by 1 + 0*218174] in the linear remainder. 
In this way we find: 

X- 1 - 0-21817 4j « 0-0508 - 0-0117J 

so that \ « 1-0508 - 0-2064J and A2 is the conju- 
gate of this. The real part of \ is very close to 
the correct value; the error in the imaginary part is 
greater, but this is not surprising since approxi- 
mations to this imaginary part are in effect obtained 
by taking the square' root of a small difference of 
relatively large quantities. 

Example 3 

Starting with the quartic equation 

(x- l)=(x- 3) ^ x^-6x3 + 12x2- iO;»; + 3 ^ (35) 

find the roots of 

X*- 6*15x3 + 12-7475x2- ll*039625x + 3*441375 

- (36) 

the exact values of which are 0*95, 1*05, 1*15 and 3. 

In the notation of equation (3): 

Sflg-^ -0*15: Sa, ^ 0*7475; 

Sai = -1*039625; Sflo ^ 0*441375 (37) 

and we then find immediately that the quantity 



called 4Kx) in equation (3) is zero for x - 3, and 
therefore SA^ - exactly. We prefer not to use any 
special technique, taking advantage of this, as our 
main object is illustration of a general method 
adaptable forprogramming; under such circumstances 
it is preferable to have as few special cases as 
possible, whereas a few extra repetitions of the 
normal routine are not particularly objectionable. 
However, the next step is to write (36) in the form: 



(x-l)3 



-0*15x=' + 0*7475x2- l*039625x + 0*441375 



X- 3 
^ +0*15x2- 0*2975x +0*147125 
- -0*000375 when x = 1 



(38) 



So the next approximations to \, Ag and A3 are 
respectively: 

X - 1 + 0*072112e(2fc + Dttj/S (k = 0,1,2) 
^ 0*927888, 1*036056 + 0*062451j 
We therefore next regard (36) as a variant of: 

{(x- 1)3 + 0*000375} {x- 3} ^ (39) 

or x'^- 6x3 + 12x2- 9.999625X + 2*998875 = 

so that Sag, Sa^ are unaltered, Sa^ - -1*04 and 
Sao ' 0*4425. Hence (36) can be written: 

(x- 1) +0*072112 

+0*15x3- 0*7475x2 + ^.q^^ - 0*4425 
^ (x- 3) {(x- 1)2- 0*072112 (x- 1) +0*005200} 

(40) 

and the next approximation to the real root is 
obtained by putting x - 0*927888 on the right hand 
side, the numerator of which in this case happens 
to be 0*15(x - 3) (x - 1) (x - 0*983333). The right 
hand side of (40) is therefore: 



0*15 X 0*072112 X 0*055445 



0*0156 



-= 0*038445 



and the corresponding approximation for a root of 
(36) is 0*9663. The corresponding equation to (40) 
when we seek to improve the approximation 1*036056 
+ 0*062451j is: 

X- 1*036056 -0*062451j WD 

0*15(x- 3) (x- 1) (x- 0*983333) 

^ (x- 3) (x- 0*927888) (x - 1*036056 + 0-062451J) 



X- 1-036056 - 0-062451J 

0-15 [(x2- 2-072112X + 1-077312) + (0-088779x- 0-093979)] 
2 X 0-062451J X (x- 0-927888) 



(42) 



When X = 1-036056 + 0-062451J and x^ - 2-072112x 
+ 1-077312 ^ 0, (42) reduces to: 

X- 1-036056 - 0-062451J ^ 0-055775 - 0-Olj 

so that X ^ 1-091831 + 0-052451], and this number 
and its conjugate are the required approximations. 
We therefore now regard (36) as a variant of: 

(x- 0-9663) (x2- 2-183662X + 1-194846) (x - 3) - 

(43) 
or 

X*- 6-149962x3 + 12-754805x2- ii.069335x 
+ 3-463739 = 

and it is now clear that we have approached the 
given equation (36) much more closely than was 
apparent by comparing the approximations found to 
the roots with the exact roots of (36). We now treat 
(36) as a variant of (43) using the procedure of 
Section 2 for well separated roots once more. We 
have, in the notation of (3): 

4ix) 

^ -0-000038x3- 0-007305x2 + 0-029710x- 0-022634 

^ -(x- 3) (0-000038x2 + 0-007419X- 0-007453) (^^j 

So the next approximation to the root near 0-9663 is 
given by: 



X- 0-9663 



0-000038x= + 0-007419X - 0-007453 



x2- 2-183662X + 1-194846 



-0-00024854 
0-018509 



^ -0-01343 



(45) 



so that X « 0-9529. The corresponding next 
approximation to the root near 1-091831 - 0-052451J 
is found, by a calculation analogous to equations 
(41) and (42), to be 1-0987 - 0-0027J. 

Now Equation (38) indicates that the initial 
approximation of three equal roots all unity is 
inevitably followed by one real root and two com- 
plex conjugates (for an equation having real co- 
efficients) whereas in fact (36) has three real roots 
(apart from x = 3). The approximation 1-0987 ± 
0-0027J obtained from (44) has a much smaller 
imaginary part than the previous approximation 
1-091813 ± 0-052451J, and this indicates that the 
next move is to regard (36) as a variant of: 



(x- 1-0987)2 (x- 0-9529) (x - 3) 

^ x"- 6-1503x3 + 12-751944x2- ll-053419x 

+ 3-450857 ^ (46) 

If we now apply the procedure of Section 3 for an 
isolated pair of equal roots, we obtain as the next 
approximation: 

(x- 1-0987)2 

0-0003x=' - 0-004444x2 + o-013794x - 0-009482 
(x - 0-9529) (x - 3) 

0-0003x^- 0-003544X + 0-003161 
X- 0-9529 

^ +0-0025422 when x ^ 1-0987 (47) 

indicating that the two roots in question are in fact 
real, and that the next approximations to them are 
1-1491 and 1-0483 respectively, which are very 
close to the true values 1-15 and 1-05. The corre- 
sponding value for the root near 0-9529 is given by: 



X- 0-9529 

0-0003x2- 0-003544X + 0-003161 




- - 


(x- 1-0987)2 

0-000056328 

— - - 0-00265 

0-02125764 


X ^ 0-9529 
(48) 



X -- 0-9663 



The next step would therefore be to regard the 
equation as a variant of: 

(x- 0-9503) (x- 1-0483) (x- 1-1491) (x - 3) 

^ x"- 6-1477x3 + 12-735891x2. ll-023107x 



+ 3-434200 - 



(49) 



and the accurate determination of the roots is now 
straightforward - the procedure of Section 2 for well 
separated roots only is required. It is at first sight 
surprising that the discrepancies between the co- 
efficients of (49) and (36) are greater than those 
between the coefficients of (46) and (36), but it is 
the values of the quantity called 4^x) in Equation 
(3) when x is equal to an approximate root of (36) 
which control the situation. 



6. ADAPTATION FOR A COMPUTER PROGRAMME 

The above examples indicate tiiat the procedure 
we have discussed should be satisfactory in all 
cases. It should also be straightforward for pro- 
gramming. The basic process involved is the 
determination of S\,- in Equation (2); this requires 
the determination of the ratio of two polynomials 
for an argument which may be complex; what is 
involved is therefore division of the polynomials by 
a linear complex divisor (or alternatively a quadratic 
divisor with real coefficients, and substitution of a 
complex number in a linear remainder). It is also 
necessary to test whether \ki - kj\ is ever small, 
and if so, how many roots are such that say: 



h- kj\< 0-1 \ki + \ 



(50) 



If there are n roots satisfying (50), then the pro- 
cedure of Sections 3 and 4 must be invoked. 

Even when this procedure does have to be 
invoked, what is needed is still the determination 
of the ratio of two polynomials for complex argument, 
and the nth root of this ratio must also be taken. 
If n - 2, the equation will not require this procedure 
more than once. If « - 3, as in Example 3, the 
special procedure may have to be invoked a second 
time for a lower value of n, usually 2, but not there- 
after. This possibility, however, is automatically 
catered for by means of the test (50). 

Subroutines for the evaluation of a polynomial 
for real or complex argument (with complex co- 



efficients if necessary), and for the determination 
of the polynomial whose zeros have known values, 
will be required repeatedly. 

The question remains of obtaining the roots of 
an initial equation which cannot be regarded as a 
variation of an equation already obtained. In the 
case of ionospheric propagation which stimulated 
the present investigation, it is known that at the 
bottom of the ionosphere, the relevant quartic 
equation has two pairs of equal roots, and therefore 
we can start calculating eigen-values at the bottom 
of the ionosphere and work upwards, even if we 
subsequently have to consider the behaviour of an 
electromagnetic wave incident at the top of the 
ionosphere and trace its progress downwards. Like- 
wise in other cases where systems of equations 
have to be considered with varying coefficients, 
there may often be one individual member of the 
system for which the equation is easily soluble by 
a special technique. 



7. CONCLUSIONS 

If a "new" algebraic equation can be regarded 
as a variation of an "original" equation whose roots 
are already known, the "new" equation can be 
solved in terms of the 'original" equation by an 
adaptation of Newton's approximation, whatever the 
relative positions of the roots of the "original" 
equation. The adaptation in question can be satis- 
factorily formulated as a computer programme. 



SMW 
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